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In this paper we study the issue of the role of nonlocality as a possible ingredient to solve long 
standing problems in the physics of black holes. To achieve this goal we analytically derive new 
black hole metrics improved by corrections from nonlocal gravity actions with an entire function 
of the order 1/2 and lower than 1/2, the latter corresponding to generalized uncertainty principle 
corrections. This lets us extend our previous findings about noncommutative geometry inspired 
black holes recently recognized as nonlocal black holes due to an entire function of order higher 
than 1/2. As a result we show that irrespective of the order of the function, nonlocality leads 
C^j ' to the following properties for black hole spacetimes: i) horizon extremization also in the neutral, 

non rotating case; ii) black hole phase transition from a Schwarzschild phase to a positive heat 
capacity cooling down phase; iii) zero temperature remnant formation at the end of the evaporation 
process; iv) negligible quantum back reaction due to the presence of an upper bound for the Hawking 
^ ' temperature. Finally we show that, in agreement with the general theory of cut off functions, a 

, regular deSitter core accounting for the energy density of virtual gravitons replaces the curvature 

Pj^ ■ singularity only in the case of entire functions of order 1/2 or higher. 
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I. INTRODUCTION 



Nonlocality is an old age theme in the quest of regular (ultraviolet complete) theories of fundamental interactions 
[l|, but for long time it has largely been ignored due to the initial success of regularization schemes. However, when 
gravity is concerned, the usual renormalization procedure runs into insurmountable difficulties and nonlocality might 
turn out to be an unavoidable requirement. This is due to the fact that we expect the spacetime to drastically 
change its nature when probed at distances comparable to the Planck length. In such a view the conventional 
| description in terms of smooth differential manifold would be just the low energy limit for a highly exited, granular 



spacetime plagued by local loss of resolution. This feature is translated as the emergence of a natural cut off, i.e., a 
" fundamental length supposed to cure the usual bad short distance behavior of gravity. Therefore it is not surprising 
that all the major formulations of quantum gravity such as string theory, loop quantum gravity, asymptotically safe 

■ gravity, noncommutative geometry, generalized uncertainty principle encode the idea of a fundamental length and 
turn out to be not local in the conventional sense. Given this background we could proceed by analyzing the nonlocal 

■ character emerging in one of the current candidate theories of quantum gravity and possibly derive phenomcnological 
f^) \ consequences from it. Even if fascinating, this route might present additional difficulties which may vary according to 

the limitations of each specific formulation. Alternatively one would like to study nonlocality, as a specific character 
by itself, irrespective of the above formulations of quantum gravity. In such a way one could hope to disclose, in a 
simpler way, deep conceptual features about the issue of the ultraviolet completion of gravity and, in principle, confront 
the conclusions with those coming from the above candidate theories of quantum gravity. This research program is 
so intriguing that many authors have employed nonlocal theories in a variety of contexts (for an incomplete list of 
examples see [2]). The standard procedure to implement nonlocal effects consists in considering an infinite order 
of derivative terms in field actions. Among the several possibilities, in this paper we want to study the effect of 
nonlocality on the gravitational field by assuming the form of the action recently proposed in |3j-l5|) 2- e -> 

Sgrav = J ^Xy/^j TZ{x) , (1) 

where in place of the Ricci scalar R(x) we have 

TZ(x) =fd 4 yV=9 A 2 (x - y)R(y) (2) 

with A 2 (x — y) a bi- local distribution defined by 

A 2 (x-y)=A 2 (n x ) S^(x-y). (3) 

Here A(O x ) is an operator of O x = £ 2 <7 M ^V A *V l/ , the dimensionless generally covariant D'Alembertian operator, and 
£ a free parameter which fixes the length scale of the theory. The above action can be seen as the leading term 
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of what recently proposed in the context of a quantum theory of gravity with higher derivative corrections to the 
Einstcin-Hilbert action |6|]. The general strategy of theories of this kind is to find a suitable form of A (D x ) so that the 
Euclidean momentum space function A(— p 2 ) plays the role of the cutoff function in the quantum gravity perturbation 
theory expanded against a fixed background spacetime at all orders. 

To reach the above goal we need to restrict the functional class where A can be chosen. For notational convenience, 
let us start by writing the Euclidean momentum space representation of A as A(z) where z — —£ 2 p 2 . According to 
Q, the first requirement for A(z) is to be an entire function when treated as a function of the complex variable z. 
Otherwise any singularities of the function A(z) for finite z will lead to the appearance of some additional non-physical 
singularities in (fTJ). Then we distinguish three possible cases. First, the function A{z) is an entire function of order 
a < 1/2, i.e., 

\A{z)\ < e°l z l° (4) 

where a is a positive number. For these functions there is not any direction in the complex plane z, along which they 
could decrease. Consequently, they cannot play the role of the cutoff functions and in this case it is impossible to 
make the perturbation theory finite. A second case occurs when the function A{z) is an entire function of order 1/2, 
i.e., 

\A{z)\ < e°l z l 1/2 . (5) 

Such a function may have only one direction in the complex plane along which it may decrease and provide a consistent 
cutoff. Finally the third choice is an entire function of order higher than 1/2, i.e., 

\A(z)\ < e ' l(|z|) (6) 

where h(\z\) is a positive function satisfying the condition h(\z\) > a\z\ x ^ 2 at \z\ — > oo for any a > 0. Such a function 
may have some regions where it converges when \z\ — > oo so that it may be chosen as a cutoff function for the 
construction of a finite unitary perturbation theory. 

In this paper we draw our attention to nonlocal gravity as emerging from virtual graviton exchange. To this purpose 
we recall that the family of noncommutative geometry inspired black holes (NCBHs) [814TTI] (for a review see [12l]) 
corresponds to black hole solutions of nonlocal gravity with an operator A represented in terms of an entire function 
of order higher than 1/2 Q, i.e., 

A(p 2 )=exp(e 2 p 2 /2). (7) 

For the strong nonlocal character implemented by entire functions of this kind, NCBHs exhibit intriguing features. 
First, NCBHs consistently reproduce at scales larger than £ the behaviour of the corresponding Einstein gravity black 
holes. Second, NCBHs enjoy regularity properties at the origin, where in place of the conventional curvature singularity 
a regular deSitter core accounts for the nonlocal, fluctuating character of the manifold in terms of quantum vacuum 
energy associated to virtual gravitons. Local violations of energy conditions at the origin confirm the departure from 
a classical description of the gravitational field at short scales. Third, NCBHs have an improved thermodynamics: 
the Hawking temperature, even for the static, neutral case, admits a maximum followed by a positive heat capacity, 
stable, terminal phase of the evaporation In addition NCBHs have companion geometries like regular collapsing 
shells [3] and traversable wormholes [l5[ . 

Despite their attractive features, NCBHs might not be the only kind of nonlocal black holes. As a consequence we 
want to study the case of entire functions different from ([7]). The goal of this investigation is twofold: on the one hand 
we would like to know if the new features of NCBHs depend on the specific choice of the above entire function or can 
be considered as general characters of nonlocality; on the other hand we would like to release the strong requirement 
of entire function of order higher than 1/2 to study milder and more economic nonlocal modifications of gravity. 
Not surprisingly, in doing so we will show that not only noncommutative geometry but also a theory based on the 
generalized uncertainty principle (GUP) is nothing else than a particular realization of a nonlocal theory. 

The paper is organized as follows: in section [TT] we derive the general form of nonlocal corrections to black hole 
metrics, in section Hill we analyze the case of an entire function of order 1/2, in section HVl we study a model of GUP 
black hole and finally in section[V]we draw the conclusions by comparing all the existing families of black holes derived 
up to now in this way. 



II. NONLOCAL MODIFICATIONS TO THE METRIC TENSOR 



To derive nonlocal effects on the gravitational interaction we need an action to describe how the nonlocal geometry, 
modeled by (JlJ, behaves when coupled to a static matter source, i.e., 

Stot = Sgrav + S'matt (8) 
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where the matter action 5 ma tt will be specified below. The field equations can be obtained by varying (|5J) with respect 
to the metric g^ v . By neglecting surface terms coming from the variation of the generally-covariant D'Alembertian, 
we find the nonlocal version of Einstein equations [3|-|5| 

A 2 (□) [R #v - \g^R^j = torGT^. (9) 

We notice that these equations can be cast in an equivalent form. The advantage of the Euclidean momentum space 
representation is that A becomes an elliptic operator and therefore has a unique, well defined inverse. One can 
conveniently work in Euclidean momentum space and rotate back to pseudo Euclidean spacetime at the end of the 
calculation. As a result one finds 

R^-\g^R = ^G% llin (10) 

where the divergence free tensor 

1^=A- 2 (D)T^. (11) 

Eq. (|10p has the familiar form of Einstein gravity coupled to a generalized source term while ([9]) describes 
nonlocal gravity equations coupled to ordinary matter T M „. The two forms are equivalent. As a static source we 
consider the case of a pressure-less static fluid at the origin whose action is given by 

d 4 x^gp (x), p {x) = -j= J (1t5(x — x(t)) (12) 

where po{x) is the energy density describing a massive, point-like particle. The corresponding energy momentum 
tensor reads 



M 

T °o = -^(r), (13) 



and the spacetime, solution of either ([9} and (fTU|) . enjoys the usual static, spherically symmetric form 

ds 2 = -f(r) dt 2 + /^(r) dr 2 + r 2 dQ 2 , (14) 

where the coefficient f(r) has the form 

/M-l-*^. (15) 
r 

The function Q(r) is the unknown and incorporates all the sought effects of nonlocality. In the limit r 3> £, the 
function Q(r) is subject to the requirement of matching Newton's constant, Q(r) — > G. Before solving the above 
equations one has to prescribe the form of A (O x ). 

III. STATIC, SPHERICALLY SYMMETRIC BLACK HOLE SOLUTION FOR a = 1/2 

We start by considering the entire functions of order a = 1/2, i.e., 

A( P 2 ) = exp(£p/2). (16) 

To determine Q{r), it is more convenient to follow the gravity field equations in the form (fTUf and derive T MZ , by 
applying A~ 2 on T^ w . For later convenience we temporarily adopt free falling Cartesian-like coordinates and we 
calculate 

T° = -MA- 2 (□) S(x) ee ~Mq{x). (17) 
As a first step we introduce the following Schwinger representation 

1 i 1 f°° 

lnA(n) = -(-ny= w ^TjJ o dss-h-<- D \ (18) 
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Second we apply the above operator on the Dirac delta, i.e. 



1 2f ^ p n 



(2tt) 3 J r \2 
As a result one obtains the action of the operator A~ 2 (□), i.e. 



A' 2 (□) S(x) = J d 3 P e-^ e l3 -P. (19) 



By integrating the r.h.s. of (fT9| one gets 



^) = ^WTW (20) 

The function g(x) is a Student's t-distribution with 3 degrees of freedom [l6j|. It admits just 2 moments. We recall 
that physically the quantity p(x) = Mg(x) is the energy density, which is smeared out for the effect of nonlocality. 
Therefore the mass M is no longer concentrated into a point but it is "diluted" according to 



M(r) = -4tt / drr 2 T° . (21) 

This mass cumulative distribution M. (r) lets us solve Einstein equations in terms of the metric component written as 
f(r) = 1 — 2A4(r)G/r. By equating the latter with (|15l) . one can determine the initial unknown 



Q(r) = 4ttG J drr 2 g(r). (22) 



The final step is to calculate the above integral. One finds 



2 , /„n 2 r/£ 

- arctan (r £) - - 

7T 7T (1 + XTl") ) 



Q(r) = G(- arctan(r/^) - - ^ , / ) . (23) 



In Fig. [TJ we have the plot of Q (r) / G. We see that nonlocal gravity matches Einstein gravity in the infrared limit 
and it is asymptotically free. This can be checked by expanding ([23l) for r ^> £ 



G(r) a G 



7r V r / 37r V r 



(24) 



which confirms that at large distances the solution (|15p tends to the Schwarzschild geometry. The above expression 
can be used to predict experimentally testable deviations to Newton's law, which turns out to be of the form 

, M ( „ f 1 mm ^ 1 



r 



r) = -G— \1+Pk< ■ ( 2 5) 



Current experimental limits on would not give tighter constraints than 1^1 /^m [17] . On the other hand we can 
expand ([23]) for r -C £ to get 



4 /r\ 3 

G(r) « - (-) G (26) 
which confirms the asymptotic freedom of the coupling. By substituting the above expression into (TT5j) . we get 



<fc 2 « - ( 1 - ^A cff r 2 ) dt 2 + + r 2 tt 2 . (27) 

This is a deSitter line element with an effective cosmological constant A e g = 8AfG£ _3 7r _1 . The regularity of the 
solution proves that our choice of entire function (|16p leads to an ultraviolet complete theory of gravity. On physical 
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FIG. 1: On the left we have the function Q{r) /G (solid line), while on the right the solid curve is the function (p + p ra d + 2px) 
(strong energy condition), the dotted curve is the fucntion (p — \p±\) (dominant energy condition) and the dashed curve is 
{p + p±) (null energy condition). In both plots we set 1 = 1. 



grounds the picture is the following. Nonlocality makes the gravitational interaction weaker and weaker as the energy 
scale increases. This is a net effect due to the emergence of repulsive quantum gravitational fluctuations in the UV 
region. The deSitter core at the origin accounts for the mean value of the energy of these fluctuations. This fact 
becomes even more evident if we study pressure terms we need to balance the generalized energy momentum tensor, 
i.e., V M *X Ml/ = 0. From the conservation equation we get 

d r %; = ~lf-\ r )d r f(r) (T/ - T °) - - r (1/ - Z°) (28) 

where %J i — diag (— p(r),p Ta d(r),p±{r),p±(r)). By solving (|28|) . we get that both the radial pressure p ra d(?*) = 
— Mg(r) and the angular pressure p±(r) = — Afg(r) — ~ M3 r g(r) turn out to be negative. As a further evidence of 
the non-classical character of the deSitter core it is useful to study energy conditions (see Fig. [1]). In a spherical cap 
around the origin the strong energy condition is violated 

IMP r 2 — f 2 

p(r)+p lad (r)+2p ± (r)= - % < 0, if r < i. (29) 

The locally exotic character of the solution is consistent with the repulsive effect of nonlocal gravity which sustains 
the generalized matter energy profile p{r) and prevents a collapse into a point like profile M6(r). 

To study the horizon equation one may wish to solve the equation /(r^j) = 0. However, for fixed £, this is 
a parametric equation depending on the mass parameter M and cannot in general be solved in a closed form. 
Alternatively, in Fig. [5J we have the function 2<f)(r)/MG (solid line), where 4>(r) = —Q{r)M/r is the new profile of the 
Newtonian potential. Horizontal lines are curves with constant ordinate —1/MG. We see that in the classical case 
(dotted line), the Newtonian potential has a hyperbolic behaviour and always admit one intersection with horizontal 
lines for any value of M. This corresponds to the event horizon of the Schwarzschild solution at r g = 2MG. Conversely 
in the nonlocal gravity case the Newtonian potential reaches an extremum (a minimum) at tq before vanishing at the 
origin. This means that intersections can occur only if horizontal lines lie above the minimum of the solid line in Fig. 
[2j More specifically there exists a value Mg for the mass parameter such that 

• for M < Mq there is no intersection, corresponding to the case of a regular manifold without horizons; 

• for M > Mq there are two intersections, corresponding to an inner horizon r_ and an outer horizon r + : 

• for M — Mq there is just one intersection, corresponding to a single degenerate horizon r^. 

By a numerical estimate we obtain Mq « 2.21 IjG and tq w 1.83 I. We notice that for large masses horizontal lines 
in Fig. [5] descend and outer horizons approach the gravitational radius r + « r g . The global structure of the solution 
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FIG. 2: The function y — 2(j>(r)/MG = —2Q{r)/(Gr) (solid line) is plotted in comparison with the classical curve y = —2/r, 
by setting 1 = 1. Intersections with horizontal curves y = —1/MG give horizon radii. 



looks like that of a Reissner-Nordstrom black hole with a regular origin in place of the curvature singularity (see Fig. 
[3]). The regularity of the manifold discloses further properties of the solution. If one consider negative values of the 
radial coordinate, one does not get a mere analytical continuation of spacetime with positive r. Since at the origin, the 
spacetime is locally fiat, the solution with positive r is geodesically complete and therefore negative r values describe 
an additional distinct spacetime. On the other hand, being arctan(r/^) an odd function, the solution with negative r 
is the mirror symmetric of the spacetime with positive r. 

To analyze the thermodynamic properties of the black hole, we perform a Wick rotation 

ds 2 = -f(r) dt 2 + r 1 (r) dr 2 + r 2 dtt 2 — ► ds% = f{r) dr 2 + / _1 (r) dr 2 + r 2 dQ?. (30) 

The periodicity of the imaginary time yields the following profile for the black hole temperature [l8| 



4vrr + V + G(r+) ) 
The entropy of the black hole can be calculated by integrating dS = dM/T and reads 



(31) 



4 J Ao G(r H ) 

where Ah — 4:irr H is the event horizon area, and A an d A + are the same quantity calculated in r and r + respectively. 
In Q31|) the term Q' /Q provides deviations with respect to the Hawking's result, which is reproduced for constant Q = G. 
By looking for the extremum of Netwonian potential <f>' — (GM/r 2 ) (1 — rQ' /Q), one obtains that the temperature 
vanishes at ro- This is in agreement with the usual thermal behavior of extremal black holes and confirms early 
expectations about a possible evaporation switching off of the Schwarzschild black hole by taking into account quantum 
back reaction effects [19| • However in our case the mechanism of horizon extremization has a different explanation: 
the switching off is driven by Q'/G which models the modifications of the nonlocal geometry at high energies. In 
such a view the Schwarzschild geometry is a "rigid" background which neglects high energy quantum modifications 
of spacetime. By calculating Q' one can display the final formula for the temperature 

T = ( 1 - 2 ; ( arctan(r + /l) - , \ ) (33) 

47T7+ ^ (l + (r + A?) 2 ) 2 V (l + (r+/£) 2 )) J K ' 

whose plot is in Fig. 2] For large r + w r g the nonlocal corrections are suppressed and one recovers the usual 
l/r+ behaviour. At r + = r max w 3.60 £, the temperature reaches a maximum T max = 1.19 x 10~ 2 t/G, before 
vanishing at r + — tq w 1.83 an extremal black hole remnant configuration. As a consequence we obtain a ratio 



FIG. 3: The solution admits one, two or no horizons depending on M. In the case of two horizons, f(r±) — 0, the Penrose 
diagram resembles the Reissner-Norstrom geometry, except for the origin where a regular deSitter core lies in place of the 
curvature singularity. 



T/M < T max /Mo « 5.38 x 10 -3 , which implies that the vacuum energy density of quantum fields propagating on the 
black hole spacetime cannot significantly modify the background geometry (fTS"|) . In other words the black hole is so 
cold that the sole term Q' /Q is sufficient to consistently describe all metric modifications during the whole evaporation 
process. 

Eq. (|32l) is a consistent generalization of the well known area law. The problem is that we cannot determine the 
exact form of the entropy analytically. However we can write (|32[) as 



1 ( A+ A Q \ 1 r+ g>(r H ) 



4 \g{r+) Giro)/ 4 J ro Q l {r H ) 
and evaluate it by an expansion in powers (ru/£) > 1 (being rjj > ro > £). As a result one finds (see (|Al|l ) 

S « ^-y A + + 0{i/r+), (35) 

which confirms our expectations: up to subleading correction, the area law is preserved in form by a mere substitution 
of the classical coupling constant G with the nonlocal one Q(r + ). 

A further inspection in the effects of nonlocality can be made by looking at the black hole heat capacity C = dMj dT, 
that can be written as 

cM-rmm-*. (36) 



rfr + J \dr 

This relation highlights the black hole exchanges of energy with the environment. Both T and dS/dr + (~ dA+/Q(r+)) 
are positive defined, while dT/dr + is negative for r + > r max , vanishing for r + = r max and positive for r + < r max . As 
a result for large r + the heat capacity tends to the classical negative value C — —2irr 2 l_/G since nonlocal effect are 
negligible at these scales. By decreasing r + the heat capacity grows negatively, departing from the usual Schwarzschild 
behaviour. At r + = r max the heat capacity has a vertical asymptote and switches sign. This means that a phase 
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FIG. 4: On the left we have the black hole temperature and on the right its heat capacity. In both plots solid lines represent 
the new behaviors while dotted lines are the classical results. In addition we set 1=1. 



transition occurs from an unstable thermodynamic phase (C < 0) to stable one (C > 0). In such a stable phase, the 
evaporation becomes a disadvantageous process, slowing down in the vicinity of the zero temperature configuration. 
From a thermodynamic viewpoint black hole remnants are asymptotic stable configurations, that cannot exchange 
energy and evaporate (T = C = 0). In Fig. ^ we see that on the left of the asymptote the heat capacity descends 
intersecting the r + -axis in r + — tq. Then the curve stops since horizons do not form at radii smaller than tq. The 
full expression for the heat capacity C(r + ) can be found in (IB1|) . 



IV. GENERALIZED UNCERTAINTY PRINCIPLE BLACK HOLES 



Generalized uncertainty relations have been postulated on the basis of string theory arguments to account the 
emergence of effective ultraviolet cutoff i.e. a minimal observable length (for a recent pedagogical review see [20(). In 
one dimension the simplest generalized uncertainty relation has the form 

AxA P > ^(l + / 3(Ap) 2 + 7 ) (37) 

where /3 and 7 are positive and independent. While in ordinary quantum mechanics Aa; can be made arbitrarily 
small, this is no longer the case for the above equation. Due to the presence of the term /3(Ap) 2 which grows faster 
than the l.h.s. term, there exist a nonzero minimal uncertainty xq in position. This result can be generalized to an 
dimensional Euclidean manifold. By postulating commutation relations like 



[q\pj] = i5)(l + ^ 2 
\pi,Pj] = 0, (38) 
one obtains that coordinate operators fulfill the following noncommutative relation 

[/?,?] = -2ip(l (39) 
where Lij is the generator of rotations. As a consequence the integration measure in momentum space 



{l + Pp 2 



\p) (p\ (40) 



becomes suppressed in the UV region [2l[ . The repercussions of GUP on the physics of black holes have been largely 
investigated (see for instance early contributions, e.g. [22I . [23|, and recent developments, e.g. [24J): the general 
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FIG. 5: On the left we have the Ricci scalar R, while on the right the solid curve is the function (p + p ra( j + 2p±) (strong 
energy condition), the dotted curve is the fucntion (p — \p±\) (dominant energy condition) and the dashed curve is (p + p±) 
(null energy condition). In both plots we set yfJ3 — 1. 



conclusion is that the introduction of the UV cut off \f]3 can prevent black holes from evaporating completely with 
a consequent remnant formation. The latter would play a significant role in cosmology as a candidate to being the 
primary source of dark matter [25| . These considerations are based on the implementation of a generalized uncertainty 
relation to particles emitted by the black holes. 

Here we want to do a step forward by deriving a black hole metric corresponding to a source term affected by the 
cut off yffi. To achieve this goal we just need to follow the above procedure: we choose a suitable entire function to 
model the suppression of higher momenta a la GUP, i.e., 

Atf) = {l + fip 2 f'\ (41) 

where the GUP cut off is (3 = I 2 . The above entire function is of order lower than 1/2, precluding the possibility of 
convergent perturbation theories. We will understand later what are the repercussions of such a low order on black 
hole metrics. As a start we need to represent the above function in terms of the D'Alembertian operator. Since our 
goal is to determine the energy momentum tensor (|11[) . it is more convenient to display the representation of A~ 2 
instead of A. By means of a Schwinger representation we write 

POO 

A- 2 (D) = / ds e -(i+(-n». (42) 



Again by employing free falling Cartesian like coordinates we can derive the action of the above operator on the static 
source, i.e., 



A- 2 (U)5{x) = /°°d S e- s ( 1+ (- v2 ))^3 Jd 3 pe^ 



1 f°° r 

/ ds e~ s / d 3 p e -*PP e iS -P 
tt) 3 Jo J 

(43) 



(2^) 3 Jo 
1 f d 3 p 



(2tt) 3 J (l + fip 2 ) 

The above integral is a well known one. It corresponds to the Fourier transform of the Yukawa potential due to the 
exchange of virtual scalars of mass ^ffi . As a result we have that the energy density reads 

M P -\ s \/-fP 

*o oB= _^- ami0 . (44) 
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FIG. 6: On the left we have the function Q(f)jG (solid line), while on the right we have the function y = 2(f>(r) / MG = 
— 2Q(r)/(Gr) (solid line) plotted in comparison with the classical curve y = — 2/r, by setting \ff3 — 1. Intersections with 
horizontal curves y = —1/MG give horizon radii. 



A first consequence of working with entire functions of order lower than 1/2 can be seen by taking the limit at short 
scales of the above profile, i.e. \x\ <C \f$- The net effect of GUP is just to replace the "classical" Dirac delta in (fT^|) 
with e - ' 31 '/^, a slightly wider distribution. However this is not enough: contrary to what happened in the previous 
case (|2T)]) . the cut off yf]3 fails to regularize the singularity in the energy density. This fact affects the structure of 
the spacetime geometry itself. By considering the conservation of the energy momentum tensor, i.e., V^T^ = 0, one 
determines the following pressure terms 

M e~ r l^ 

2W(r) =-j-^r (45) 

^ = -2?— + 8^ e ^ (46) 

with — di&g(—p(r),p ra _d(r),p±(r),p±(r)) as above. This permits to evaluate the curvature scalar R — SirGI^ 
by tracing Einstein equations (|10p 

The curvature singularity is not removed despite the introduction of the cut off \ffi and the usual local violation of 
energy conditions (see Fig. (|5]l). In any case the GUP effects modify the profile of the singularity a fact that discloses 
some new features. By integrating (1441) one determines the profile of the cumulative mass distribution M.(r) as well 
as the running coupling Q{r), i.e., 

M(r)/M = G(r)/G = 1 - er r '^ - (r/V^)e" r/v/? . (48) 

In Fig. ©, we see that the coupling constant becomes weaker at short scales even if not enough to improve the 
manifold. Correspondingly A4(r) is not enough spread out to cure pathologies. In any case at the level of Newtonian 
potential, the usual —l/r singularity is tamed: by considering the expansion at short scales of cf> = —GAi(r)/r (or 
equivalently —Q{r)M/r) we find 

cj>{r) » —GM (?-), r<^fp (49) 



while GUP corrections quickly die off at scales r ^> y^8. This behavior of the Newtonian potential has important 
consequences on the thermodynamic side. We recall that zeros of the Hawking temperature coincide with stationary 
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FIG. 7: On the left we have the black hole temperature and on the right its heat capacity. In both plots solid lines represent 
the new behaviors while dotted lines are the classical results. In addition we set yffi = 1. 



points of the Newtonian potential. Falling off to zero at both r — > and r — > oo, the function <j) must admit a 
maximum a fact that permits to conclude that the black hole will not evaporate off but will rather condensate into 
a zero temperature remnant. This can be viewed also in Fig. ([6]): at ro 1. 73-^/3 there exists a lower bound 
Mo « 1.68V/3/G for the mass spectrum below which no horizons form. We conclude that the global structure of the 
resulting solution 



ds 2 



(l - GAf 7 (2;r/^)/r) dt 2 - (l - GM 7 (2; r/y/p)/r) 



dr 2 



l dn 2 



(50) 



is analogue of the Reissner-Nordstrom geometry with the number of horizons determined by the value of the mass 
parameter M with respect to Mq. Here j(2;r/y/p), the incomplete gamma function defined by 



x dt 



t a e 



(51) 



provides an equivalent representation of the function Q(r)/G in f|48|) . We notice that, even in the absence of charge, 
the solution develops a naked singularity for M < M . 

By following the procedure of the previous section we now display the black hole temperature 



T = 



1 



1 - 



9 

r± e 



P 7(2; r/y/jS) 



(52) 



Being T vanishing r + — > oo as well as r + — > ro it must admit a maximum, as is evident in Fig. ([7]). Qualitatively the 
temperature has a behavior in agreement to what found in (22, 25] for r + > r max s» 4.20t/J3 with r max = T(r max ) ~ 
1.89 x 1Q~ 2 ^J]3/G. Our result is a little bit more sophisticated: by recalling that stationary points in the temperature 
correspond to asymptotes in the heat capacity we obtain a phase transition to from a Schwarzschild phase (C < 0) 
to a cooling down phase (C > 0) with final evaporation switching off (C = T = 0). As in the case of the previous 
section quantum back reaction effects turn out to be negligible. This can be seen by comparing the energy of the 
emitted particles, i.e. T, with the black hole mass M: the ratio T/M < T max /Mo ~ 1.13 x 10~ 2 is so small that the 
spacetime (|50|) remains unmodified during all the evaporation process. 

To confirm the above reasoning we need to determine the black hole heat capacity. As a preliminary step we display 
the black hole entropy by using the general equation (|32l) 



S 



2tt 
~G 



dr 



H 



7(2; th/VP)' 



(53) 



Even if we cannot integrate the above formula analytically, it is not hard to prove that at the leading order for 
(rn/y/P) > 1 the entropy fulfills the area law, i.e., S ~ 7rr^/G7(2; rjyfft). On the ground of ([55)1 we can use (|3"o| to 
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determine the heat capacity that reads 

c= G7( w^) T(g(r+) - 4 ' rrr '- (54) 

where the function H(r + ) is defined by 

7(V+/vp) 7(2;r+/v/5) V7( 2 ;r+/vp)/ 

with the derivatives of the gamma functions 

7 '(2;r + /v^) = je- r+/VP (56) 
7 " (2; r + /v^) = I e -WV? (l - ^) . (57) 

The term (H(r + ) — AnT) — 4nr + (dT/dr + ) determines the sign of C according to our expectations. In particular 
when (dT/dr + ) = (i.e. at r — r max ) the heat capacity has a vertical asymptote and the black hole undergoes a 
phase transition. The resulting plot is given in Fig. [7] 



V. CONCLUSIONS 



In this paper we have analyzed the possibility for nonlocality to be the key ingredient for curing the conventional 
pathologies of classical gravity at short scales. To reach this goal, we have derived black hole solutions by extending 
our studies on nonlocal gravitational actions with entire function of order a > 1/2 0, 0, [Io| to the case of entire 
functions of order a < 1/2. Specifically we have derived a black hole solution corresponding to the case a = 1/2 and 
another one corresponding to GUP corrections modeled in terms of an entire function belonging to the class a < 1/2. 
We have now a more complete scenario and we are ready to summarize our results by pinpointing the general features 
of these black holes. 





ro[<| 


M [£/G] 


''max 


T m ax[^/G] 


(T/M) < 


Newton <f> 


curvature 


thermodynamics 


"NC-Schw" 


3.02 


1.90 


4.76 


0.0150 


7.89 x 10~ 3 


regular 


regular 


locally stable 


"NC-Dirty" 


3.02 


1.90 


4.84 


0.0148 


7.79 x 10~ 3 


regular 


regular 


locally stable 


"a = 1/2" 


1.83 


2.21 


3.60 


0.0119 


5.38 x 10-" 


regular 


regular 


locally stable 


"GUP" 


1.73 


1.68 


4.20 


0.0189 


11.25 x 10-* 


regular 


singular 


locally stable 



TABLE I: Summary of parameters and properties of some families of short scale improved black hole metrics. Specifically we 
have the noncommutative Schwarzschild black hole ("NC-Schw") the noncommutative dirty black hole ("NC-Dirty") 

- both corresponding to nonlocal gravity solutions with entire function of order higher than 1/2 - the nonlocal gravity black 
hole with entire function of order 1/2 ("a = 1/2") and the GUP modified black hole ("GUP") - corresponding to a nonlocal 
gravity solution with entire function of order lower than 1/2. 

As expected from the general theory Q , nonlocality determines a UV complete theory of gravity only if the order 
of entire functions is a > 1/2. We showed that this can be translated into the "language of black holes" by saying 
that curvature singularities are replaced by regular deSitter cores which account for the quantum vacuum energy 
of virtual gravitons. In the case of a < 1/2 such a vacuum energy develops (a fact that can be viewed by looking 
at the violation of energy conditions in the vicinity of the origin) but it is not regularized. The net effect is just 
a deformation of the profile of the curvature singularity whose degree of divergence is decreased. As a result for 
all a one finds that the Newtonian potential admits a maximum and vanishes at the origin. On the other hand at 
large scales nonlocal effects quickly die off and all black hole metrics consistently match the classical Schwarzschild 
geometry. We showed that the presence of stationary points of the Newtonian potential for all a has repercussions in 
two directions: i) the global structure of the corresponding black hole solutions is modified by allowing the horizon 
extremization; ii) the Hawking temperature does not diverge because the black hole undergoes a phase transition 
from a negative heat capacity evaporation phase to a positive heat capacity cooling down which terminates with a 
zero temperature remnant configuration. Accordingly the black hole temperature never reaches values high enough to 
disrupt the geometry by permitting a consistent semiclassical description of all the evaporation process. The remnant 
mass is not only the mass of the black hole in the extremal configuration but also the threshold mass for developing 
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horizons: lighter masses lead to horizonless regular geometries for a > 1/2 and naked singularities for a < 1/2. In 
addition for all a the area law of the black hole entropy is preserved in form: the only modification occurs at the level 
of coupling constant where the usual G is replaced by the nonlocal running coupling Q (r + ) . 

In summary we showed that the way nonlocality deforms the short distance behavior of gravity has model indepen- 
dent characters with the only exception of the singularity avoidance in the case a < 1/2. This fact opens the route 
to further considerations and new investigations. Specifically we can extend some recent results obtained for the case 
a > 1/2 to the case of generic a. Concerning the phenomenology of microscopic black holes in particle detectors 
nonlocal deformations of black hole thermodynamics might lead to distinctive signatures. Even if in the present paper 
we studied the four dimensional case only we argue that the possibility of having a maximum black hole temperature 
would imply a Hawking emission in terms of soft particle on the brane, a fact that is in marked contrast to the results 
for Schwarzschild black holes as shown in (2(| [27j . On the other hand in the presence of a non-asymptotically flat 
term, i.e. a cosmological constant, nonlocal black holes have a thermodynamics equivalent to a Van der Waals gas in 
which the role of the average volume occupied by a particle is played by the cut off I (or equivalently v7? m the GUP 
case) [ll|. On the ground of this result, one may think to exploit the black hole phase transition to describe, via the 
gauge-gravity duality paradigm, the nuclear matter transition from the "classical" skyrme crystal to a phase where 
quantum effects weaken the bond between baryons and eventually melt the crystal into a nuclear gas with consequent 
switch of the deconfinement phase into a cross over (for more details see [28| )■ As a related issue to black holes we 
recall that the UV completeness of these theories can lead to a two-dimensional Planck scale spectral dimension for a 
diffusive process propagating over a nonlocal background manifold [29j . This fact is important since it lets us conclude 
that nonlocality might be the key ingredient to disclose self-renormalizability properties of gravity at the Planck scale, 
by virtue of a mechanism of dynamical dimensional reduction. 

All the above considerations hold at a qualitative level while the families of nonlocal black holes are quantitatively 
slightly different (see Tab. Ufor a summary of their properties and parameters). We think that their differences are 
an advantage rather than a limitation. For instance having different values for the threshold mass implies a variety 
of scenarios for the production of black holes, that, by considering only the case of a > 1/2 theories, would be out of 
reach of the present ground based experiments |3Cf . In an ideal world one would like to reverse the logic by trying 
to determine the exact profile of the entire function A from experiments. This can be done also in the absence of 
cxtradimension: it has been shown that neutrino oscillations are sensible to nonlocal gravity modifications opening 
the possibility of direct observations in neutrino telescopes like IceCube and ANTARES (3lj . 

We conclude by saying that nonlocal deformations of gravity offer reliable solutions to long standing problems in 
quantum gravity and a rich phenomenology with distinctive predictions that might be potentially observable in a near 
future. 
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Appendix A: Useful formulas 

The expansion of the entropy (|35[) in powers of x = rjijl > 1 can be obtained by the following relation 
G'(x) . , s , , f , 9 s-l 1 1 \ / , . 1 



-A H {x) = tt 2 x 2 2 (x 2 + 1) - 2 + 4 — —2 arctan (x) + 



Q 2 {x) " w V x 2 (l + x-2)' x *(l + x -*f)\ x(l + x~ 2 ), 

« 16.t- 2 -32z- 4 + — -L +48x- 6 - — +0 (x~ 8 ) . (Al) 

3 7T X 15 7T X y ' 
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Appendix B: Black hole heat capacity 



The heat capacity in (|36|) reads 



i 



47r(r+/f) 



-kQ{t. 

2G 



-6( r+ /ey 



A(r + /£f 



(l + {r + /t?) (l + (r + /if) (l + {r + /lf) 



(Bl) 



2G 



1 



4 7 r(r+/£) i 



(47rr+T) 
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